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Abstract. We study in this paper an M/M/l queue whose server rate de- 
pends upon the state of an independent Ornstein-Uhlenbeck diffusion process 
(X(t)) so that its value at time t is fj,<f>(X(t)), where <j>(x) is some bounded 
function and > 0. We first establish the differential system for the condi- 
tional probability density functions of the couple (L(t), X(t)) in the stationary 
regime, where L(t) is the number of customers in the system at time t. By 
assuming that (f>(x) is defined by o(z) = 1 — e((x A a/e) V (— b/e)) for some 
positive real numbers a, b and e, we show that the above differential system 
has a unique solution under some condition on a and b. We then show that this 
solution is close, in some appropriate sense, to the solution to the differential 
system obtained when is replaced with = I — ex for sufficiently small e. 

We finally perform a perturbation analysis of this latter solution for small e. 
This allows us to check at the first order the validity of the so-called reduced 
service rate approximation, stating that everything happens as if the server 
rate were constant and equal to fi(l — eS,(X(t))). 



1. Introduction 

We consider in this paper an M/M/l queue with a server rate varying in time. 
We specifically assume that the server rate at time t is equal to /j,(f>(X(t)) for some 
function </> and some auxiliary process (X(t)). Throughout this paper, we shall 
assume that the modulating process (X(t)) is an Ornstein-Uhlenbeck process with 
mean m > 0, drift coefficient a > 0, and diffusion coefficient a > 0. This process 
satisfies Ito's stochastic equation 

(f) dX{t) = -a(X(t) -m)dt + adB(t), 

where (B(t)) is a standard Brownian motion. The stationary distribution of the 
process (X(t)) is a normal distribution with mean m and variance a 2 /(2a); the 
associated probability density function is defined on the whole of R and is given by 

, , , > dof. 1 fa ( a(x 

(2) n(x) — — \ — exp 

a V 7T 



Throughout this paper, we shall assume that the Ornstein-Uhlenbeck process is 
stationary. 

If L(t) = j denotes the number of customers in the M/M/l queue and X(t) = x 
at time t, then the transitions of the process (L(t)) are given by 

J j + 1 with rate A, 
3 ~* 1 j - 1 with rate fx4>(X(t)). 



Key words and phrases. M/M/l queue, Self-Adjoint Operators, Perturbation Analysis, Power 
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In the following, we shall assume that the condition p =' X//i < K(<f>(X (0))) < 1 is 
satisfied so that it is straightforward to show the existence of a stationary proba- 
bility distribution for the Markov process (X(t), L(t)); see Meyn and Tweedie [T2] 
for example. 

The study of the above system is motivated by the problem of bandwidth sharing 
in telecommunication networks and the coexistence on the same transmission links 
of elastic traffic, which adapts to the level of congestion of the network by achieving 
a fair sharing of the available bandwidth, and unresponsive traffic, which consumes 
bandwidth without taking care of other traffic. See for instance for a discussion 
about bandwidth sharing in packet networks. The choice of an Ornstein-Uhlenbeck 
process as modulating process is natural for several reasons: mathematically this is 
a standard "typical" diffusion process with an equilibrium distribution and secondly 
it can be seen as a centered approximation of the number of jobs of an M/M / oo 
queue (the unresponsive traffic), see for example Borovkov [4] or Iglehart [6] or 
Chapter 6 of Robert [17] . 

One of the objectives of this paper is to investigate the so-called Reduced Service 
Rate (RSR) property for which the system would behave as if the server rate were 
equal to the mean value /jX.((/)(Xq)) . Even though some results can be established 
for arbitrary perturbation functions <j){x), we shall pay special attention in the 
following to the case when the function 4>{x) has the form 

(3) cj>(x) = 1 - s((x A (o/e)) V (-b/e)) 

for some small < e < 1 and real numbers < a < 1 and b > 0, where we use 
the notation nVi = max(a, b) and a A b = min(a, b). The choice of the bounded 
perturbation function is discussed at the end of the paper. 

As it will be seen, one of the important technical problems encountered in the 
perturbation analysis is the existence of a reasonably smooth density probability 
function for the couple (X(t), L(t)) in the stationary regime. Conditions on e for 
ensuring the existence and the uniqueness of a density probability function will be 
established in the following via Hilbertian analysis. More precisely, let Pj(x) denote 
the stationary probability density function that the process (L(t)) is in state j and 
the process (X(t)) is in state x and let P denote the vector whose jth component 
is Pj(x)/ p 3 . In a fist step, we show that P is solution to an equation of the type 

(4) Qf + V(4>)f = Q, 

where fl is a selfadjoint second order differential operator on 2?' (R) N , D'(M) denoting 
the set of distributions in K . Unfortunately, the operator V((f>) is not selfadjoint 
so that Kato's perturbation theory for selfadjoint operators cannot be applied. 
Nevertheless, we prove that the above equation has a unique non null smooth 
solution P E C 2 (]R) N for sufficiently small e. In addition, we prove that when 
replacing 4>(x) with = 1 — ex, we obtain an equation of the type 

(5) n/ + eV($)/ = 0, 

which has a unique non null smooth solution for sufficiently small e, V{$>) being 
independent of e. Denoting this solution by g, we prove that P and g are close 
to each other for some adequate norm when e is small. We then perform a power 
series expansion in e of g and we determine the radius of convergence of this series. 
By explicitly computing the two first terms of the series, this eventually enables us 
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to prove the validity of the reduced service rate approximation at the first order for 
the system. 

The problem considered in this paper falls into the framework of queueing sys- 
tems with time varying server rate, which have been studied in the queueing liter- 
ature in many different situations. In Nunez-Queija and Boxma |15j . the authors 
consider a queueing system where priority is given to some flows driven by Markov 
Modulated Poisson Processes (MMPP) with finite state spaces and the low priority 
flows share the remaining server capacity according to the processor sharing disci- 
pline. By assuming that arrivals are Poisson and service times are exponentially 
distributed, the authors solve the system via a matrix analysis. Similar models 
have been investigated in Nunez-Queija [131 114) by still using the quasi-birth and 
death process associated with the system and a matrix analysis. The integration of 
elastic and streaming flows has been studied by Delcoigne et al. [5], where stochastic 
bounds for the mean number of active flows have been established. More recently, 
priority queueing systems with fast dynamics, which can be described by means of 
quasi birth and death processes, have been studied via a perturbation analysis of 
a Markov chain by Altman et al pQ. A probabilistic analysis of these queues with 
varying service rate has been presented in Antunes et al. [21 [3] . Our point of view 
in this paper is completely different since a functional analysis approach is used to 
tackle the perturbation analysis. The key difficulty for the case considered in the 
present paper is that the associated Markov chain has an infinite state space. 

This paper is organized as follows: In Section [21 we establish the basic system of 
partial differential equations for the joint probability density functions of the pro- 
cess (X(t), L{t)). We recall in Section [3] some basic results on the generators of the 
Ornstein-Uhlenbeck process and the occupation process in an M/M/l queue. In 
Section[4]it is proved that this system has a unique solution with convenient regular- 
ity properties in an adequate Hilbert space when e is sufficiently small. In Section [5l 
wo carry out a perturbation analyoio for the perturbation function <f»(x) = 1 — ess-, 
we show that when replacing <j>(x) with <j>(x) = 1 — ex, the corresponding differ- 
ential system has also a smooth solution in the underlying Hilbert space when e 
is sufficiently small. We then prove that the solutions to the differential systems 
for <fi defined by Equation ([3]) and <E> are close to each other in some appropriate 
sense. By expanding the solution of the second differential system in power series 
of e, we show that at the first order the so-called Reduced Service Rate property 
for the original system holds; the subsequent terms of the associated expansion are 
also expressed. Some concluding remarks are presented in Section [5J 



2. Fundamental differential problem 

2.1. Notation and differential system. The goal of this section is to establish 
the fundamental differential system for the conditional probability density functions 
Pj{x), j > 0, in the stationary regime, where Pj(x) is the probability that the process 
(L(t)) is in state j knowing that the Ornstein-Uhlenbeck process (X(t)) is in state 
x. As long as we do not have proved regularity results, these functions have to be 
considered in the sense of distributions, i.e., for all j > 0, Pj(x) G 2?'(R), where 
V(M.) is the set of distributions in M. The distributions Pj[x), j > 0, are formally 
defined as follows: for every infinitely differentiable function with compact support 
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(p(x) (denoted, for short, ip € Cg°(R)) 

Pj(x)n(x)<p(x)dx = - / ip , (x)¥(X(0)<x,L(0)=j)dx, 



where n(x) is the normal distribution given by Equation ([2]). 

Throughout this paper, we shall use the following notation. The functional space 



L 2 (R, n) = <^ / : R -> R : / f{x) 2 n{x)dx < oo 

I JR 

is a Hilbert space equipped with the scalar product defined for /, g 6 L 2 (K, n) by 

(i ',3)2= / f{x)g{x)n{x)dx 



and the norm of an element / € L 2 (R,n) is H/H2 = v/ (/, 3)2- If W is a separable 
Hilbert space equipped with the scalar product (., .)^ and associated norm ||.||h> 
we define the Hilbert space 



L 2 (R,i l ;H) = |(/ J (i),j>0)a 2 (M;ii)™: / \\f(x)\\f i n(x)dx < 00 
equipped with the scalar product 

if ,9)= I {f{x),g{x))-Hn{x)dx 



Finally, let I 2 (p) be the Hilbert space composed of those sequences (cj , j > 0) taking 
values in R and such that Y^jLo ° 2 jP^ < 00 ' anc ^ equipped with the scalar product 
defined by: if c = (cj) and d — (dj) in l 2 {p), (c, d) p = X^^Lo c j < ^jP^^ ^ ne associated 
norm is defined by: for c e £ 2 (p), \\c\\ p = yj (c, c) p . Let ej denote the sequence with 
all entries equal to except the jth one equal to 1. The family (ej,j > 0) is a basis 
for I 2 (p) . The space £\ (p) denotes the subspace of I 2 (p) spanned by the vectors ej 
for j > 1. 

In a first step, we determine the infinitesimal generator of the Markov process 
(X(t), L(t)) taking values in 1 x N. We specifically have the following result. 

Lemma 1. The process (X(t), L(t)) is a Markov process mKxN with infinitesimal 
generator Q defined by 



a 2 d 2 f df 
(6) Qf(x,j) = Y^2"( x ^') ~ a ( x ~ m )~dz( x >ti 

+ X (f(x,j + 1) - f(x,j)) + p<j>(x)l {j>0} (f(x,j - 1) - f(x,j)) , 

for every function f(x,j) from R x N in R, twice differentiable with respect to the 
first variable. 

Proof. According to Equation JT]), the infinitesimal generator of an Ornstein-Uh- 
lenbeck process applied to some twice differentiable function g on R is given by 

(?) Hg= ^^ [x) _ a{x _ m) ^ [x) . 

The second part of Equation ^ corresponds to the infinitesimal generator of the 
number of customers in a classical M/M/l queue with arrival rate A and service 
rate p(fi(x), when the Ornstein-Uhlenbeck process is in state x. □ 
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For finite sequences / = (fj{x)) of infinitely differentiable functions with compact 
support, the equation for the invariant measure for the Markov process (L(t), X(t)) 
is given by J2j>a IrS f(x, j)pj(x)n(x)dx = 0, that is, 

^ f°° (a 2 d 2 f 1 , .dfc 
f^ J-oo V 2 dx 2 dx 

+/J,l{ j>0 }<f>(x)fj-i(x) - (A + fj,(f>(x)l {j>0} )fj(x) + Xf j+ i{x) ) pj-(a;)n(a;)da: = 0. 

Via integration by parts, we obtain for every finite sequence / = (/j(x)) of infinitely 
differentiable functions with compact support 



a 2 d 2 Pj dP 3 



4-i \./_oo \ 2 dx 2 dx 

j>0 



a(x — m )—r^ ) fj(x)n(x)dx 



DC 



(nMj>o}<t>(%)fi-i{x) ~ {^ + ^(t>{x)t{j > o})fj{x) + \f j+1 {x))p : j(x)n(x)dx 

= 

and then 



3>0 



r°° fgicfPj dPj 
7-oo V 2 dx 2 dx 



+ (j,l {j>0 yPj-i(x) - (A + ( u<?!>(x)l{ i> o})Pj(x) + \<f>(x)Pj + i(x) ) fj(x)n(x)dx = 0. 
This implies the following result. 

Proposition 1. 77ie family (Pj(x) == pj{x)/ p* ,j > 0) G 2?'(]R) N is solution in the 
sense of distributions to the following infinite differential system: for j > 0, 

+ til {j>0 }Pj-i(x) - (X + ^(x^^Pjix) + Xi>(x)P j+1 (x) = 0. 

2.2. Additional properties. For the system considered in this paper, we have 
p4>(x) > — a) for all a; € K and < a < 1. Classical stochastic ordering 
arguments imply that the process (L(i)) is stochastically dominated for the strong 
ordering sense by the queuing process of the M/M/l queue with input rate A and 
service rate — a). Hence, if the solution (Pj(x),j > 0) of the infinite differential 
system (|SJ) is related to the conditional probability density functions (pj(x),j > 0) 
of the couple (X(0), L(0)) as Pj(x) = Pj{x)/p> for all j > 0, then 

(9) Va;eR,Vj>0, Pj(x) < 

If a < 1 — yfp, it is easily checked that for all i£l, the sequence (Pj(x),j > 0) is 
in the Hilbert space £ 2 (p). 
In addition, for all j > 



(10) / pj(xyn(x)dx < P(L(0) = j) < oo, 

Jm 

since pj(x) = P(i(0) = j \ X(0) = x) < 1. It follows that for all j > 0, the function 
Pj (x) should be in the space L 2 (R, n) . 
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Hence, if the solution (Pj(x),j > 0) of the infinite differential system is 
related to the conditional probability density functions (pj(x),j > 0) of the cou- 
ple (X(0), L(0)) as specified above, then Pj(x) £ L 2 (M.,n) for all j > 0. From 
inequality (jTUJ) , we also deduce that if a < 1 — y/p, 

since F(L(0) = j) < pP /{I - a) j . 

If follows from the above remarks that to show the regularity of the conditional 
probability density functions pj (x) for j > under the assumption a < 1 — y/p, we 
are led to prove that the differential systems admits a unique regular solution in 
the space L 2 (K, n; l 2 {p))- 

In the next section, we review some properties of the operators associated with 
the generators of the Ornstein-Uhlenbeck process and the Markov process describing 
the number of customers in an M/M/l queue. 



3. Some results on the operators associated with the generators of 
the Ornstein-Uhlenbeck process and the M/M/l queue 

It is well known in the literature (see for instance [TH]) that the operator H 
defined by Equation ([7]) is selfadjoint in the Hilbert space L 2 (R, n). The eigenvalues 
of this operator are the numbers —aj, j > 0, and the normalized eigenvector 
associated with the eigenvalue —aj is the function hj given by 



(11) 



hj(x) 



1 



H j (V2a(x - m)/cr), 



where Hj(x) is the jth. Hermitc polynomial. The sequence (hj,j > 0) is an or- 
thonormal basis of L 2 (M.,n). The domain of the operator H is the set 

D{H) = {/ E H 2 {R,n) : x 2 f G L 2 (R,n)} , 

where H 2 (R,n) is the Sobolev space defined as follows: 

H 2 (R, n) = {/ € C^R) :/,/'€ i 2 (M, n) and the weak derivative /" e i 2 (R, n)}. 

While the operator H io well known in the literature, lcao information is available 
on t The operator A associated with the Markov process describing the number of 
customers in an M/M/l queue and defined in £ 2 (p) by the infinite matrix 



(12) 



.4 



/-A A 







A 


(i 


-(AH 







V • 







A 

-(A + /i) A 



has already been studied in the technical literature, notably in [7] (see also [5]). 
From these references, we know that the operator A is selfadjoint. Associated with 
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the operator A is the operator A\ defined in £ 2 (p) by the infinite matrix given by 



/-(a + m) 



(13) 



A* 





A 

-(a + m) 

A* 





A 

-(a + aO 

A* 





A 
-(A + /x) A 



V 



Note that the above matrix is the generator of the Markov process describing the 
number of customers in an M/M/l queue and absorbed at state (see [7]). Finally, 
let A^ denote the truncated operator associated with the finite matrix 



(14) 



A[ 



[N] _ 



/-(a + m) 

A 1 




V 



A 

"(A + A*) 





A 

-(a + aO 





The above matrix is the generator of the Markov process describing the number of 
customers in the finite capacity M/M/l/N queue and absorbed at state 0. (See 
[5j for a related model.) In the following, we recall the spectral properties of the 
operators A, A\ and A^; see [7j|8j|9] for details 

Lemma 2. The operator A in t 2 (p) is bounded and symmetric, and then selfadjoint. 
In particular, for all f £ £ 2 (p), 



(15) 



o<(-A/,/) p < m (i + vp) 2 II/II 



which implies that the operator —A is monotonic (i.e., (—A/,/) > for all f £ 
£ 2 (p)). There exists a unique normalized measure dip(z), referred to as spectral 
measure, whose support is the spectrum o~(A) of operator A, and a family of spaces 
{H z (p)}, z £ cr(A), such that 

• the Hilbert space £ 2 (p) is equal to the direct sum of the spaces 7i z (p), i.e., 
every f £ £ 2 (p) can be decomposed into a family (f z ,z £ a (A)), where 
f z £ 1~lz(p) and J \\fz\\pdip(z) < oo. Moreover, 



(f,9) f 



(fz,9z) P dip(z) 



• The operator A is such (Af) z — zf z for z £ o~(A), where (Af) z is the 
projection of (Af) on the space H z (p). 

The spectral measure dij){x) is specifically given by 



(16) / h(x)dip(x) = (1 - P )h(0) 

VP 




2VV 



dx, 



- M (l + V/3) 2 

for any smooth function h. 

The spectrum of the operator A is a {A) = [— p(l + ^/p) 2 , — p{\ — y/p) 2 ]U{0}. The 
operator A has a unique eigenvalue equal to 0, and the eigenspace H-o(p) is spanned 
by the vector e with all components equal to 1. For z £ (— (vA+ \fp) 2 ^ ("\/A — 
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^/p) 2 ), the space 7i z {p) is the vector space spanned by the vector Q(z), whose com- 
ponents Qj{z), j > are defined by the following recursion: 



(17) 



Qo(*) = l, Qi(z) = + A)/A 

fiQ j+1 (z) - (z + A + p)Qj{z) + M<3j-i( z ) = 0, j > 1. 



T7ie vectors (Q(z)) for z G (— (VA + ^/m) 2 , — (^/m— %/A) 2 ) /orm an orthogonal family 
with weight function dip(z): for all j, k, 

Qj(x)Q k (x)dip(x) = -^S jtk , 
pj 

where Sj^k is the Kronecker symbol, equal to 1 if j = k and if j =/= k. 

Note the polynomials Qj{x) appearing in the above result are known as per- 
turbed Chebyshev polynomials in the literature on orthogonal polynomials [18) . 
For the operator A\, we have the following result, where we use Chebyshev poly- 
nomials of the second kind (U n (x)) defined by the recursion 



(18) 



U (x) = l, Qi(x) = 2x 

U j+ i(x) = 2xUj(x) - Uj-i(x), j > 1. 



Lemma 3. The operator A\ in the subspace span(ej, j > 1) is bounded and sym- 
metric, and then selfadjoint. In particular, for all f G span(ej, j > 1), 

(19) M (i - Vp) 2 II/IIp < Mi/, f) P < m(i + VpfWfWl- 

The associated normalized spectral measure dipi (z) is given by 



(20) #r(x) = ~Wl ( X+ 2 ^ ) l { , e( -Mi+yp) 2 ,-,(i-VP) 2 )} ^ 



The spectrum of the operator A\ is diffuse (there are no eigenvalues) and equal to 
the interval a (Ai) = [— p(l + a//?) 2 , — /i(l — \fp) 2 ] ■ The Hilbert space sp&n(ej,j > 1) 
is equal to the direct sum of the spaces 'Hz(p), for z G (— p(l + ^/p) 2 , — p(l — ^fp) 2 ), 
where the space H z (p) is the vector space spanned by the vector Q^\z), whose 
components Q^\z), j > are defined by the following recursion: 



(21) 



Q ( o\z) = l, Q{ 1 \z)^(z + \ + p)/\ 

pQfl.iz) - (z + A + n)Qf\z) + nQf\{z) = 0, j > 1. 



The vectors (Q^(z)) for z G (— (vA + ^fji) 2 ,— {^fpZ — \/A) 2 ) form an orthogo- 
nal family with weight function dipi(z). The polynomials (Qj(z)) are related to 
Chebyshev polynomials as follows: for j > 0, 



Q^(z)-—U- ( Z + X + V ) 



Finally, for the operator A^\ we have the following result. 
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Lemma 4. The operator is symmetric (and then selfadjoint) in the vector 

space span(ei, . . . , ejv) equipped with the scalar product induced by (., .) p . The eigen- 
values of the operator A± are the solutions to the polynomial equation 



1 (-x;l) = Q$(-x;l), 



where the polynomials Qj(x) are defined by the recursion (|2ip . The eigenvalues 



are denoted by 



, N with x± 



[N] 



<4 JV] 



< 



< X 



[N] 
N 



The vec- 



tors Q^' N \xj), j = 1,... , N , form an orthogonal basis of span(ei, . . . , ejq), where 
Q( 1 ' N )( X j) is the vector with the kth component equal to Q^l (xj), k, j = 1, . . . , N, 
where the polynomials (Q^\z)) are defined by Equation (|21[) . 

The operator A naturally induces in L 2 (R, n; l 2 (p)) an operator that we still 
denote by A. The same property is valid for the operators A\ and A^ in the 
spaces L 2 (R, n; span(ej, j > 1)) and L 2 (R, n; span(ej ■, j = 1, . . . ,N)), respectively. 
Similarly, the operator H induces in L 2 (R, n; l 2 (p)) an operator that we still de- 
note by H and which is defined as follows: for / 6 L 2 (R,n;l 2 (p)), Hf is the 
element with the jth component equal to Hfj. This operator also induces in 
L 2 (K, n; span(ej, j > 1)) and L 2 (R, n; span(ej, j — 1, . . . , iV)) operators denoted by 

Hi and h\ , respectively The operators A, Ai, A^, H, Hi and are clearly 
selfadjoint in the spaces where they are defined. 

With the above definitions, the fundamental differential system 



(22) (H + A)f + Vf = 0, 

where the operator V is defined by: for / 6 L 2 (R, n; I 2 (/?)), Vf 
where B is the operator associated with the infinite matrix 

■ ■ \ 



reads 



(<Mx) - l)Bf, 



B 



In the notation of Equation ((4]) , we have f2 
(-b/e))B. 



fo 


A 








-A* 


A 








-M 











V 








A 







= H 



The matrices 



Bi 



and 



B[ N] = 





A 










-fJL 


A 


. 










A 











—fi A 


(-, 


A 












A 















V • 










A and V(e) 
\ 



-e((x A (o/e)) V 







define the operators Bi and in the spaces L 2 (R, n; span(ej ;, j > 1)) and 



L (R, n;span(ej,j = 1 

the operator £?, i?i and _B| JV| are bounded with a norm less than or equal to 



,N)), respectively. It is straightforwardly checked that 
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+ ^/p)- With the operators Bi and Bi are associated the operators V\ and 

induced by V in the spaces L 2 (R, n; span(ej, j > 1)) and L 2 (R, n; span(e,, ■, j = 
1,...,N)), respectively. 

In the next section, we prove that the differential system (J8J) (or equivalently 
Equation (|22|) ) has a unique solution in L 2 (R, n; I 2 (p)) . 

4. Existence and uniqueness of a solution 

When we refer to the existence of a density probability density function satisfying 
the differential system ((8|), we think of a vector (Pj(x)) such that every function 
Pj{x) is twice continuously differentiable over R (i.e., Pj(x) G C7 2 (R) for all j > 0). 
But, this differential system may have a solution, which is in L 2 (R, n; £ 2 (p)) but 
with components not in C 2 (R). In the following, we prove that the differential 
system (|8]) has a unique solution in L 2 (R, n; i 2 (p)) and then we show at the end of 
the section that the components of the solution are C 2 (R) functions. 

If / e L 2 (R,n;( 2 (p)) is solution to Equation ((22]), then 

(23) (H 1 +A 1 )f 1 + Vif 1 = -fie 1 {f ), 

where f 1 is the projection of / on the space L 2 (R, n, span(ej ■, j > 1)) and ei(/o) 
is the element of L 2 (R,n, span(ej, j > 1)) with the first component equal to fo(x) 
and all other components equal to 0. 

The operator (Hi + A\) is self-adjoint and invertible. It is straightforwardly 
checked that the norm 

\\(H 1 +A 1 )- 1 \\ = f 

supllKi/i + Ax)- 1 /!! :/ e£ 2 (R,n,span(e,,j> 1)),||/|| =1} 

is such that 

since by Lemma [3] and the monotonicity of the operator —Hi, we have for all / in 
L 2 (R, n, span(ej, j > 1)) 

(-(Ht + A 1 )f, f) > (-H x f, f) + (-Ax/, /) > -mi - Vp) 2 II/II 2 - 

In addition, if / G L 2 (R, n; i 2 (p)) is a solution to Equation (f22|) . then the function 
S j=o P 3 fi ( x ) e L 2 (n) since by Schwarz inequality 

2 

1 POO 00 -j^ 

n(x)dx < — — / \fji x )\ 2 P jn ( x ) dx = T— ll/ll 2 < °°- 

*■ P J -CO j_Q J- P 

By summing all the lines of Equation (|22[) . we see that the function X)jlo P 3 fi( x ) 
has to be solution to the equation Hg = in L 2 (R, n). 

By using the above observations, we prove the existence and uniqueness of a non 
trivial solution in L 2 (R, n; £ 2 (p)) to Equation (f2"2")l by showing the following results: 
(1) If e satisfies some condition, Equation (f2"3"]) has a solution for all fo £ L 2 (n); 
if (fj) is the solution, we set 

oo 

(24) Kf Q = J2p i fj&L 2 (R,n). 
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(2) If / is a non trivial solution to Equation (j22|) . then /o is a non trivial 
solution to equation 

(25) f + Kf = l. 

(3) Equation (f2"5)) has a unique non trivial solution. 

To prove the last point, we intend to use the Fredholm alternative, which requires 
that the operator K is compact. 

Lemma 5. Under the condition 

(26) aVb< v , 

where a V 6 = max(a, 6), Equation (|23[) ftas a unique solution in L 2 (IR, n; £ 2 (p)). 

Proof. Equation (|23[) can be rewritten as 

(I + (Hi + A 1 )~ 1 V 1 )f = -fi(Hi + A 1 )- 1 e 1 (f ). 

Since V x is bounded so that for all / g L 2 (R, n; l\(p)), |(Vi/,/)| </i(aVi))(l + 
\/p)II/II 2 ' we deduce that the operator (Hi + A 1 )^ 1 V\ is bounded with a norm less 
than or equal to (a V b)(l + v / p)/(l — -*//)) 2 • Under Condition (f2"6"l) , the norm of 
operator (i/i + Ai)" 1 ^ is less than 1 and we then deduce that (I+(H\ + -Ai) _1 Vi) 
is invertible [16] and Equation (|23|) has a unique solution in L 2 (R, n; t\(p)). □ 

It is worth noting that under Condition (|26|) . we have a < (1 — v/p)- Moreover, 
the above result ensures that the operator K is well defined by Equation ([24]) . Now, 
we prove that under the same condition, the operator K is monotonic. 

Lemma 6. Under Condition (|26[) . the operator K is monotonic, which implies that 
there exists at most one non trivial solution to Equation ()25[) . 



Proof. We first note that from Equation (|24[) . we have 

(A7o,/o) 2 - (-fi(H 1 +A 1 + V 1 )- 1 e 1 (f ),e(f )), 

where e(/o) is the vector with all entries equal to /o. The above equation can be 
rewritten as 

(Kf , foh = (-m(I + V(#i + Vei(/o), e(f )). 

Since —pA^ei(fa) = e(/o), we have 

(27) / ) a = ((I + V(tfi + U 1 ))- 1 e(/ ), e(/ )) 

and hence (Kf , f ) 2 > 0. Indeed, for all / g L 2 (M, n; ^f(p)) 

((i + Ar^i + ^vo/,/) = ((i + Vffi^/R^rVr/,/) 

/ (avb)(i + yp) \ 
- v, (i-Vp) 2 / 11/11 ' 

where we have used the fact that the operator A~^ l H\ is monotonic, ||Vi|| < p(a V 
b)(l + y/p), and H^i" 1 !! ^ l/(A t (l — V/') 2 )- ^he a bove inequality implies that 

(Ot+V^i+V^) -1 /,/) 

> fi- (a l m ±f> ) II (i + + ^VO-Vf > o, 



(i-VP) 2 

and Inequality ([27]) follows. □ 
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We now turn to the compactness of the operator K. The major difficulty comes 
from the fact that the operator (Hi + Ai)^ 1 is not compact, since we know that 
the spectrum of this self-adjoint operator is not discrete. However, by truncating 
the infinite matrix defined by Equation (|12|) , we can introduce compact operators 
and subsequently prove that the operator K is compact. Let us fix some N > 0. 
We first prove the following technical lemma. 

Lemma 7. For N > 0, the operator + A^)^ 1 in L 2 (R, n\ span(ei, . . . ,ej\r)) 

is compact. 

Proof. By using Lemma [4] and the orthonormal basis (h n ) is L 2 (R, n), we know 
that the family ej.k{x) = hk(x)Q^ 1 ' N ' (xj) for k > and j — 1,.,.,N forms an 
orthogonal basis of L 2 (R, n; span(ei, . . . , ejv))- In particular, we have 



(H[ N] +A[N] 1 )- 1 e J Ax) = - 



ka - 



(x). 



The operator (H[ 

operators 
l,...,iV,fc = 0, 



[JV] 



) then appears as the norm limit as M — * oo of the finite 

M) by 



rank operators (h[ N ' M ^ + A^' M ^) 1 defined in the vector space span(ej fc(ic), j 



(h[ n < m] +a[ n ' m] ) 



(x) 



,[N] 



ej,k(x) 



and the result follows. 



□ 



Lemma 8. Under Condition J26J), the operator K is compact. 

Proof. Let us consider a bounded sequence (/q) in L 2 (M, n). (Without loss of 



generality, we assume that H/0II2 = L) Since the operator (H[ ' + A[ ') 1 V{ 1 is 
bounded with a norm less than or equal to (a V b)/j,(l + J~p)/xi , and the operator 

(H[ N] + A 1 ^)- 1 is compact, we deduce that the operator (h[ N] + A l f ] + V^)' 1 
is compact. 

Let /* = (fj(x)) denote the vector — ii(H\+ Ai + Vi)^ 1 ei(f^) . Since the operator 
(H\ + A\ + is bounded with a norm less than or equal to l/(/x((l — ^fp) 2 — 

aVk(l + \f~P))): the vector (/*) is such that 



(28) 

We have 



N 1 
\pt U 1 - VP) - fflV6 ( 1 + v^)) 




"M/o 




In view of inequality 



the sequence appearing in the right hand side of the 



above equation is bounded. Hence, since the operator (h[ N] + A[ N] + V^y 1 is 
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compact, it is possible to extract a sub-sequence (/ l 



[N] 







such that 



V -A^ + i + J 



( /r 

J2 

J* CO 

V In 



as fc — > oo in L 2 (K, n; span(ei, 
hand side of the above equation is in L 

fik . fx 

>N JN 



, ejv)), where the vector appearing in the right 



(R, n; span(ei, 
In particular, we have — * in L 2 (R, n) as A; 



■ - ,ejv)). 

oo. This implies that 



/ /iv 
f' k 



N+l 

■ik 

N+2 



I 



[JV] 



[JV] 



A 



[N] 







DC 

JV 



as fc — > oo in L 2 (R, n, £ 2 (p)), since the operator (i?-^ 
We set 

/ /at+i 



A [JV] + v , 1 ,v J) _ 1 . g boundei 




r[JV] 



OO 

N+2 



= { h[ n] + a[ n] + v} n] )- 1 



\ 



The vector with the jth component equal to f°° is in L 2 (M.,n,£ 2 (p)) and we have 
Kfo° ~ ¥ Yl'jLiP'fj m L 2 (R,n) as k — > oo. We then deduce that from every 
bounded sequence (/q) in L 2 (R, n), we can extract a sub-sequence (/q*) such that 
K/q 1 is converging in L 2 (R, n). The operator X is hence compact. 



□ 



By using the above lemmas and the Fredholm alternative, we can state the 
following result. 



Proposition 2. Under Condition (|26[) . Equation (|25p has a unique solution. This 
establishes that Equation (|22[) has a unique non trivial solution in L 2 (R, n, £ 2 (p)) 
if Condition (|26j) is satisfied. 



The above result has been established for the perturbation function cj)(x) = 
1 — e((x A (a/e)) V (— (b/e))) and we have exploited the fact that the function 
|1 — 4>{x)\ is bounded by a V 6. In fact, it is possible to prove a similar result when 
4> is replaced with Q(x) = 1 — ex. Let us define the operator W in L 2 (R, n; £ 2 (p)) 
by: if /=(/, (*), j>0) 

Wf = -exBf. 
Note that the domain of W is given by 



D(W) = { (f 3 (x),j>0)eL 2 (R,ny 



/ x 2 fj(x) 2 n(x)dx < oo > 

3=0 J 



In the notation of Equation ([5]), we have V($>) — —xB. 

We can then state the following result, whose proof is given in Appendix |A"1 
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Proposition 3. Under the condition 
the equation 

(30) (H + A + W)f = 
has a unique solution in L 2 (K, n, £ 2 (p)). 

To prove the existence and the uniqueness of the solutions to Equations l]22p 
and (j3"0|) , we have only supposed that the components of the solutions are in D(H), 
in particular the components are in H 2 (R, n). But, by examining the differential 
systems satisfied by the different components, notably by taking into account the 
continuity of the functions <p(x) and $(x), these components are clearly in C 2 (R). 

To conclude this section, let us mention that the solution to the differential 
system ([8]) have the following probabilistic interpretation: for all j > 0, in the 
stationary regime Pj(x) = ¥(L(t) — j \ X(t) — x)/ pi . In particular, we have for all 

i>o 

(31) 0<P,(z)<l 
and for all x G M 

oo 

(32) Y^l'r >■>!>■ I. 

j=0 



5. Perturbation analysis 

The goal of this section is to prove the following Reduced Service Rate approxi- 
mation. 

Theorem 1. For sufficiently small e, the first order expansion of the generating 
function of the stationary distribution of (L(t)) is given by 



E 



1 B rr^-Frf m£+ " (£) - 



Therefore, K(u L ^) ~ E(w Le ), where L e has the stationary distribution of the 
number of customers in an M/M/l queue when the server rate is 1 — em. This 
shows a principle of reduced service rate approximation, i.e., everything happens 
as if the server rate were fixed equal to 1 — me. 

To show the above result, we proceed as follows: 

(1) We compare the solutions to Equations (|2"2"]) and (|30[) when Conditions ([2^]) 
and (|29[) are satisfied. In particular, we compute an upper bound for the 
norm of their difference. 

(2) We develop the solution g to Equation (|30| in power series expansion of e. 
In particular, we explicitly compute the two first terms. 

(3) We finally prove Theorem [T] 

Throughout this section, we denote by P and g the solutions to Equations J22| 
and {30]) in L 2 (M, n; t 2 (p)), respectively. 
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5.1. Comparison of the solutions P and g. The solutions to Equations |22|) 
and (|30|) when they exist are close to each other when e is small. We specifically 
have the following result. 



Proposition 4. Assume that Conditions (|26[) and (j29|) are satisfied. These solu- 
tions P and g are such that 

(33) \\P-g\\<D(e), 

where the function D(e) is given by 



1.3-1 1 Di:) = M[\ + x Iy^~ p + X Jy^~ p M ) A( :) 



with 



A(e) 2 = (fj 2 + 3\n) ((ex-a) 2 l {x > a/e} + (ex + b) 2 l {x <_ b/£} )n(x)dx. 

JR 

The function D(e) is O (e 5/2 exp(-a(a A fe) 2 /(2cr 2 e 2 ))) when e ^ 0. 
Proof. The vector P is such that 

(H + W + A)P + (V - W)P = 
and hence, if P 1 denotes the projection of P on the space span(ej,j > 1), 
(Hi + Ai + W^P 1 = (Wi - VijP 1 - Mei(Po) 

and then 

(35) P 1 = {H 1 +A 1 +W 1 )- 1 (W 1 -V 1 )P 1 -fi(H 1 +A 1 + W 1 )- 1 e 1 (P ). 
It follows that 

oo 

Y^PiW = (P\e 1 ) p = ((H 1 +A 1 +W 1 )- 1 (W 1 -V 1 )P 1 > e 1 ) p 

3=1 

- M ((Pu + A x + W 1 )- 1 e 1 (P ), e 1 ),. 

We have 

+ A! + V^O-^^Po), e^p = X' 5o , 

where the operator K' is defined as the operator K (defined by Equation (|24|) ) but 
by replacing V with W. In addition, since P satisfies (P, e) p = 1, we come up with 
the conclusion that Pq verifies 

1 - P = ((Pi + Ai + Vi)- 1 ^! - yi^.e'jp + iC'P . 

Since go + X '90 = 1, we obtain 

(P - ,g ) + K'(P - g Q ) = ((£Ti + Ai + Wi)" 1 ^ - Vi)PV)p- 

Since the operator fsT' is monotonic, the above equation implies that 



IP0-.90II2 < jj^—um + A^w^n^-v^p^ 
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The norm \\(Wi — V\)P 1 \\ is given by 
\\(W 1 -V 1 )P 1 \\ 2 = 

oo „ 

V / ((ex - aft {x > a/e} + (ex + b) 2 l {x <_ b/e} ) (nPj(x) - XP j+1 (x)) 2 p j n(x)dn 

and by using Equation pip , we obtain 
\\(W 1 ~V 1 )P 1 \\ 2 

<(p, 2 +3\p) / ((ex - a) 2 l {x > a/e} + (ex + b) 2 l {x <^ b / e} ) n(x)da 

JR 

Simple computations show that 

\2 n 



5e 5 cr 5 



{ex -a) l {x > a/s} n{x)dx ~ —^j 

5e 5 cr 5 

{ex + b) 2 l {x <^ b/e} n(x)dx 



-e 



8a 5 / 2 b 3 ^ 

whene -> 0. The term \\{Wi~V 1 )P 1 \\ is hence <3(exp(e 5 / 2 exp(-a(aA&) 2 /(2cr 2 e 2 ))) 
when e — > 0. In addition, since by using Equation (|6ip in Appendix [Al 



M (i - vp) 2 V (i - Vp) 2 

we deduce that ||/o — go\\2 is dominated by a term, which is 0(exp(e 5 / 2 exp(— a(aA 
b) 2 /(2a 2 e 2 ))) when e ->■ 0. Finally, by using the fact that (H 1 + A x + Wi)g = 
—pei(go), we deduce from Equation ([55)) that 

P 1 - g 1 = (H x + Ai + W 1 )- 1 (Wi - Vi)P 1 - KHi + A 1 + W 1 y 1 e 1 (P - g ), 

which implies that 

WP'-g'W < \\{Hi + A 1 + W 1 )- 1 \\\\(W 1 -V 1 )P 1 \\ + X\\(H 1 +A 1 + W 1 )- 1 \\\\P -g \\2 

and hence IjP 1 — g 1 \\ is dominated by a term which is 0(exp(e 5 / 2 exp(— a(a A 
b) 2 /(2a 2 e 2 ))) when e -> 0. □ 

5.2. Power series expansion of the solution g. 

5.2 A. Notation. We assume that the solution g to Equation ([30[) can be uniquely 
decomposed as a power series expansion of the form 

(36) g = gW+eg^+e 2 g^ + ...., 

where S L 2 {R, n; £ 2 {p)) for z > 0. In addition, to facilitate the computations, 
we shall consider the generating function g u {x) = YlJLo 9j( x ) u:> P 3 \ which is an 
element of L 2 (M., n; £ 2 {l/p)). Indeed, if g is written in the form 

oo 

where Cj(u) — J2m=o c j-m um with (cj tm ,m > 0) € (- 2 (p), then g M (a:) can be written 
as 

oo 

9u{x) = ^2C J {u)h j {x), 
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with Cj(u) = J2m=oCj,mU m = cj(pu). Since (cj, m ,TO > 0) € £ 2 (p), (Cj. m ,m > 
0) e t 2 {l/ P ). Finally, we have ||g|| 2 = Nl? = ££o HQIIi/ P - 

The generating function g u (x) will be expanded as 

(37) 9u (x) = g(0) ix)+£g il) {x)+£ 2 g (2) {x) + . . . . 

where € L 2 (]R, n) for all i > 0. The function gu\x) corresponds to the case 
e = and is given by 

(0)M - C (g) 



(38) gW(x) 



1 — pu 



where c(g) is the normalizing constant. 

In the following, we prove that the elements g^ have to satisfy a recurrence 
relation of the form g^ = Q^xg^ -1 ^) for i > 1 and for some linear operator 
whose norm is finite. 

In the following, we assume that the expansion (|3T[) is valid and we investigate 
the conditions which have to be satisfied by the elements gW. In a first step, we 
prove the following property satisfied by the functions (gu\x)). 

Lemma 9. For i > 0, the vector j s j n L 2 (M, n; £ 2 (p)), where Lf{R, n; i 2 (p)) is 
the sub-space o/L 2 (R, n; £ 2 (p) composed of those elements (fj(x)) such that fj(x) £ 
span(/io, . . . , hi) for all j > 0, i/ie functions hj being defined by Equation (jlip . - </ie 
function gffl (x) in Expansion (|37[) hence satisfies for N > i 

(39) lim -L fl W(a.)=0. 

x — >±oo x 

Proof. The proof is by mathematical induction. The result is true for i = since 
g(o) _ c (^) e ( e being the vector with all components equal to 1). 
If the result is true for i. From Equation (|30[) . we have 



(H + A)g {l+1) = -Wg {i) . 
By using the recurrence relation satisfied by Hermite polynomials |10j 
(40) H j+1 {x) - 2xH j {x) + 2iH j _ 1 (x) = 0, 

it is easily checked that the image by the multiplication by x of span(/io, . . . ,hi) is 
span(/i , . . . ,h i+ i). Therefore, since by assumption g( l > belongs to L 2 (R,n;£ 2 (p)), 
we immediately deduce from the uniqueness of the decomposition on the basis 
(hiej,i > 0,j > 0) of the Hilbert space L 2 (K, n; £ 2 (p)) and the selfadjointness of 
the operator H + A, that is in L 2 +1 (R, n; i 2 (p)) and the result follows. □ 



5.2.2. First order term. In a first step, we pay special attention to the derivation 
of the first order term because it gives the basic arguments to derive higher order 
terms. Moreover, the explicit form of the first order term will be used to examine 
the validity of the reduced service rate approximation (see Theorem [l| . 

On the basis of the domination property given by Lemma [9l we explicitly com- 
pute the function (x). From Equation (f30|) . it is easily checked that the function 
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9u\%) satisfies the equation 



-*(H( rf>w+ ^) 



where the constant v(u) is given by 

- u){l - pu) 

v{u) = . 

au 

In a first step, we search for a particular solution to the ordinary differential 
equation 

c 2 <9 2 £ u d£ u p(i(l-u)c(g) 

of the form 

£«(ai) = a(u) + b(u)x. 
Straightforward manipulations show that 

a{v{u) - l)c[g) v(u) 

Noting that ^oO^) = 0, it follows that if we write gu(x) = £ u (x) + ip u (x), then 
the function ipu(x) is solution to the equation 

(42) ~2~dx^ ~ a< " X ~ m ^~tx + ay ^ u )^ u i x ) = P (~ ~ 1 ^) ^o(^)- 

By using the domination property of Lemma El we can determine the form of the 
function ipo( x )- 

Lemma 10. The function tpo(x) is given by 

y/a(x-m) 

Wo{x) = c + ci 

a 

for some constants cq and C\ . 

Proof. By introducing the function k u (x) defined by 

, v ( ctix — to) 2 \ , . 

(43) k u (x) = exp f \ a2 J fa (x) 

and then the change of variable 

(44) z = ^ x - m \ 

a 

Equation (|4"2")l becomes 

(45) ^ + {2u{u) + 1 - z 2 )k u {z) = *H(l-i) k (z). 

oz z a \u J 

The homogeneous equation reads 

d 2 k u 



dz 2 



+ (2u(u) + 1 - z z )k u = 0, 
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which solutions are parabolic cylinder functions (see Lebedev [10] for details). Two 
independent solutions vi(u; z) and v 2 {u: z) of this homogeneous equation are given 
in terms of Hermite functions as 

2 In 2 try 

(46) vi(u;z) = e~ z 1 H v ( u )(z) and v 2 {u; z) = e z 1 H_ u ( u )_i(iz). 
The Wronskian W of these two functions is given by 



W(z) 



- „-(v+l)7Ti/2 



By using the method of variation of parameters, the solution to Equation l|45p 
is given by 

feu (■*) = li{u)vi(u\ z) + ~f 2 (u)v 2 (u; z) 

_ (I _ l") eKDWJ f [ Vl ( u; y) V2 ( u - z ) _ Vl (u- z ) V2 (u- y) ] ko (y) dy : 
a \u J J 

where 71(11) and 72 (u) are constants, which depend upon u. 

The function ijj u {x) enjoys the same domination property as function g^p(x), 
given by Lemma [9] Hence, for TV > 1 



(47) 



lim -^e* l 2 k u {z) = 0. 

z— >±oo z 



From Lebedev [TO], we have the following asymptotic estimates 



(48) 



H v {z) ~ (2 Z ) 1 



.fc=0 



it! 



when \z\ — > 00 and | argz| < 37r/4 — 6 for some 5 > 0. Moreover, when z — > —00 



r(-f) 



ur"-v 



n 1 

£-( i , + l) 2fc (2z)- 2fc + 0(|zr 2 "- 2 ) 



.fe=0 



N 



(2zf , k£N. 

The above asymptotic estimates and Lemma © imply that for u g (0, 1) such that 
v(u) £ N with v > 1, we have 

71 («) = -— - e ( ' +1)7 " /2 r v 2 {u;y)k {y)dy 
= -^(--l) e^W 2 v 2 (u; y)ko(v) dy 



a \u 



and 



72 (u) 



a \ ii 



The latter equation implies that for all n > 1 



(49) 



- y2 / 2 fc (y)-ff42/)rfy = o, 



where H n (x) is the nth Hermite polynomial. By Property (|47jl . the function 
y — > exp(y 2 /2)fco(y) is in L 2 (R, exp(— y 2 ) dy). Since Hermite polynomials form 
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an orthogonal basis in this Hilbert space, Equation (|49|) entails that the function 
y — ► cxp(y 2 /2)fco(y) is orthogonal to all Hcrmitc polynomials H n with n > 1 and 
then that this function belongs to the vector space spanned by Hq and H±. Hence, 
function ko(z) should be of the form 

k Q (z) = (c + ciz) e - z2/2 
for some constants cq and c± and the result follows. □ 

By using the above lemma, we are now able to establish the expression of g2~ \x). 

Proposition 5. The function g-u (x) is given by 

(50) gl 1) (x) = 



(1) _ c(g)(u! - - pu)(u- 1) 



u\U\(l - p)(u - ui)(l - pu\)(l - pu) 2 

c(g)(l-u) 



(u - Si)(l - pUi)(l - pu) 
where u\ and U\ are the two real solutions to the quadratic equation 



pu 2 - !•//•-) - I = 



(V 



with < u\ < 1 < U\ 

Proof. By taking into account Lemma [TU1 the function K u {z) defined by 
and the change of variable ([44]) , satisfies the equation 



(51) ^ + (2u(u) + l -z 2 )K u {z) 



a \ u J \ 1 — pu \ y/a 



We search for a particular solution of the form 

K u {z) = {a{u) + b(u)z) e~ z2/2 . 

Straightforward computations yield 

/ \ 1 f . Pucjg) 
a{u) = — r c + m 



(1 — pu) \ 1 — pu 
b{u) = a - ») _ ( Cl , P™<9) 



p(u - ux){u - ui) \ Va(l -pu) / 
It follows that the general solution to the above equation can be written as 

(52) K u (z) = (a(u) + b(u)z) e~ z2/2 + 7i(u)«i(u; z) + j 2 (u)v 2 (u; z), 

where the functions v\ and v 2 are defined by Equation (|46|) and the constants 71 (u) 
and 72 (u) depend upon u. 

By differentiating once Equation (f5"2")l with respect to z and using the fact that 
the Wronskian of the functions Vi(u; z) and V2(u; z) is exp[(i>(u) + l)7ri/2], we can 
easily express 71 (u) and 72(11) by means of K u (z), a(u), and b(u). This shows that 
7i(u) and 72 (u) are analytic in the open unit disk deprived of the points and u\. 
From the asymptotic properties satisfied by the functions v\ and v 2 , we know that 
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71 (it) = and 72 (u) = for u such that v(u) > 1. It follows that 71 (u) = 72(1*) = 
for \u\ < 1. 

By using the fact that gi 1 ^ (x) has to be analytic in variable u in the unit disk, 
we necessarily have 

pauic(g) 
^/a(l - pui) 

and then, 

ac(g)(l - u) 



b(u) 



Moreover, since g[ 1 \x) = 0, we have 



y/a(u - Gi)(l - - pit) ' 

pc(g)m 



Co : . 

1-p 

and then, 

p(u - l)c(g) 



(53) a(w) - (l-p)(l-pu)3-- 
By using the expressions of a(u) and b(u), the result follows. □ 

5.2.3. Higher order terms. We assume that gi l \x) can be expressed as 

i 

(54) 9 { i ) {p)=Y J ^)hj{^ 

where the function hj is defined by Equation (fTT|) and the coefficients Cjj are 
analytic functions in variable u. This assumption will be justified a posteriori. 
From previous sections, this representation is valid for z = 0,1. If it is valid for 
i — 1, then the function g„ (x), i > 1, satisfies the equation 

(55) — — — a(x — to)— \- av{u)g y > {x) 



2 <9x 2 9a; 



First note that by using the recurrence relation ([40]) satisfied by Hermite poly- 
nomials, it is easily checked that 

i 

x(gt 1] (x) - gt 1] (x)) ='£d itj (u)h j (x), 

3=0 

where 

di i(u) = ^— =(cj_i - Cj_i ,i_i(0)), 
2va 

and for < j < i — 1, 

= 7T7^( c i-i,j-i( w ) - c i-i,i-i(0)) + mfe-i^it) - Cj-x^O)) 

+ -^fe-ij+iW ~ Ci_i J+ i(0)). 
By using the above notation, we have the following result. 
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Proposition 6. The coefficients Cij appearing in the representation (|54p of gu\x) 
are recursively defined as follows: we have 

i \ 1 ~ P 
1 — pu 

and for i > 1, 

di,o(u) - di, (l) 



Ci,o( u ) = 



1 — pu 



where for j > 1, Uj and Uj are i/ie two rea/ solutions to the quadratic equation 
v(u) = j, i.e. 

pu 2 - (l + p + u + 1 = 

with < Uj < 1 < Uj . 

Proof. As in the previous section, we first search for a solution to the equation 



n 2 «2<rOQ atM 



Assuming that the function ffl (x) is of the form 



k=0 

we have, by using the fact that the functions hj(x) are eigenfunctions of the oper- 
ator H associated with the eigenvalues —aj and that these functions are linearly 
independent, for j = 0, . . . , i, 

a \u J v(u) - j 
It is easily checked that £q (x) = 0. We can then decompose g$ (x) as as 

where the function ipffl (x) is solution to the equation 



2 dx 2 



(x- m)^- + v{u)^\x) = p (~ - lj ^(x). 



By using the same arguments as in the proof of Lemma [TU1 we can easily show that 
%1)q\x) has the form 

i 

3=0 
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where the coefficients Cj e C for j = 0, It follows that the function (x) is 
solution to the ordinary differential equation 

Y~£hk~ ~ a ( x ~ m )~£r + ai/ ( u )9u } ( x ) = M [ ~ - 1 ) z2( c J + d i,j( u )) h j( x )- 

By using the same arguments as in the proof of Proposition [5l we come up with 
the conclusion that gu\x) is of the form (|54|) with the coefficients Cjj-(u) given by 

a \u J v{u) — j 

Since the function <?! (x) has to be analytic in the open unit disk, we have for j ' > 1 

cj = —dij(uj) 

In addition, since g\ (x) = 0, we have cq — —dij(l). □ 
The normalizing constant c{g) is chosen such that 

gi(x)n(x)dx = 1. 

3 

From the above analysis, we see that gf 1 (x) = for i > 1 so that c(g) = 1 — p. 

5.2.4. Radius of convergence. In this section, we examine under which conditions 
the expansion ([36]) defines an element of L 2 (M.,n,£ 2 (p)). In a first step, note that 
as a consequence of Proposition [6l the function gu (x) can be written as 



where the operator 8 is defined in L 2 (R, n; £ 2 (l/p)) as follows: for an element 
/ e L 2 (R, n;£ 2 (l/p)) represented as 

oo 

Mx) d = y^Cj{u)hj{x), 



3=0 



the element F — 9/ is defined by 



DC 



^— ' V u J v(u) — 7 

where we set Uo — 1 and u = 1/p. 

It is easily checked that for j > 1, < Uj < 1 < < Moreover, the 

function Cj(tt) appearing in the expression of /„ is analytic in the disk D p = {z : 
\z\ < 1/y/p} and continuous in the closed disk D p = {z : \z\ < 1/y/p} for j > 0. 
Similarly, for all j > 0, the function 



- j 

is analytic in D p and continuous in D p . With the above notation, we can state the 
main result of this section. 

Proposition 7. The operator O is bounded and if e < l/(m||0||), where ||0|| 
denotes the norm of O, then the sequence defined by Equation (|36[) ( or equivalently 
by Equation (|37j) ) is in L 2 (R,n;£ 2 (p)). 
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Proof. Let / e L 2 (R, n; £ 2 (l/p)) be defined by the function 

oo 

fu{x) = ^C {u)h 3 {x). 
3=0 

For (cj) e l 2 (l/p) associated with the generating function 

oo 

C ( U ) = J2 C 3 U ^ 
3=0 



we have 



1 



c i /„ = 



1 



d0. 



Let us moreover define the sequence (cj ) associated with the generating function 

c(u) - 

Assume first that j > 1, then 



c(u) = -(1 — u) 



p / 1 _ \ c(m) - c(Mj) 
a\« / K u ) - j 

1 c(tt) - c(uj) 



and then 



VpJ (uj - i/vp) 2 27r Jo 



2tt 



e i6 lsfp~-Uj 



d6. 



Simple manipulations show that 



1 

2^ 



2tt 



o i0 



IVp- 



d9 < lie" 2 



1 + 



1 



1 — yO?i 



It follows that ||c||i/ p < Kj||c||i/ P , where 



- 1 + ^ 



VpJ _ 1 /Vp)( 1 /Vp- u j) 
i + V~p 



1-pu 2 



(i - vp) 2 + 



i 



l - pu 2 



It is easily checked that the sequence (kj) for n > 1 is decreasing. 
When j = 0, we define 



c(u)-c(l) _ c(u)-c(l) 



1 — /Olt 

It is then easily checked that ||c||i/ p < «o|| c lli/pi where 

1 



K = 



Vp 



1 + 



1 

~ ~p 



Define k = maxj^o, Ki}. The above computations show that for all / £ £ 2 (l/p), 
II < k ll/H • It follows that the operator O is bounded; its norm is denoted by 
||9|| = inf{c > : V/ e H, ||0/|| < c||/||}. The above computations show that 



(56) 



|6|| 



1 + 



1-p 
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From the above computations, we deduce that 

H fl (0|| < ||ein| c (°)"|| 

where the sequence c^* 1 is associated with the function 

x . 

1 — pu 

Straightforward computations show that 

( o).i„2 = H2 . ( ^ 



(57) 11^11=^; h - 

where Hi(x) is the ith Hermite polynomial. Using the asymptotic estimate 
we have 

||c (0) *i ~ m l 

when i — ► oo . It follows that ||cW|| < a, with a, ~ (||6||m) 1 as i tends to infin- 
ity. It follows that the sequence defined by the expansion ([571) is convergent in 
L 2 (R,n;£ 2 (p)) if e||6||m < 1. □ 



5.3. Proof of Theorem[H We assume that Conditions (|26|) and ([29]) are satisfied. 
By observing that E (u L W) = (P, U), we have for u < 1 



E 



- (g, U) < \\P - g\\ \\U\\ = -^—\\P g\\ < -^—D(e) 
V / 1 — pu 1 — pu 



where U is the vector of L 2 (R, n, i 2 {p)) with the jth component equal to u J and 
D(e) is defined by Equation (13~3| . 

From the power series expansion of g, we have 

oo 1 oo 

(g, U) - (g(°\U) - e(g^,U)\ < ]T \\U\\ < ]T «CjV , 

with 

where we have used Equations ([55|) and (|Ff|) . We clearly have 



1-P 

f — pu 

(9 {1 \U) = |°° g^(x)n(x)dx = a(u) = - 

where a(u) is defined by Equation ([53")) . Theorem [1] then follows. 

To complete the analysis, note that the generating function of L S) the stationary 
number of customers in an M/M/l with input rate A and service rate fi(l — me), 
is for \u\ < 1 and when e < (1 — p)/m 

(o) m _ r„(D rrv = p(i~">" 2 £ 2 



E{u L ')-(g^,U)-(g^,U) 



(1 - p - me)(l - p)'< 
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Hence, by gathering these relations and by taking u = 1, we obtain an uniform 
bound for the difference between E (u L ^) and E (u Le ) for u £ [0, 1]. 



(58) sup 

0<u<l 



< E B 



dcf. 



1-P 



D(e) 



1 



1 / \ i 2pm 2 

— > k(j)£ j + 7--, rr, 

yj ' (i _ p - me)(l - pf 



J=2 



Below are some numerical experiences on the role of e and a on the bound Eb for 
a = 1/2, b = 1 rn = 1, X = 7, a = 1 and p = 10. It is reasonably low for small 
values of e, it seems to be quite sensitive on the values of the parameter a as the 
figures below show. 




(7 = 2 

0.0002 0.0004 0.0006 0.0008 " 0.001 ° 

FIGURE 1. The bound e — > E B of Relation (JSHJ) for er = 2, 3, 4. 



x = 5 




x = 3 



x = 2 



x = 1 



50 a 



FIGURE 2. The bound a — * of Relation ((SSJ) for e = x • 10 
with x = 1, 2, 3, 5. 



PERTURBATION ANALYSIS OF AN M/M/l QUEUE 



27 



6. Concluding remarks 

The perturbation analysis performed in this paper has allowed us to prove the 
validity of the so-called reduced service rate approximation for the system consid- 
ered under some specific conditions. Such an approximation is very important from 
a practical point of view because each type of traffic can be considered in isolation, 
the impact of unresponsive traffic on elastic traffic is only via the mean value. 

The results presented in this paper have been obtained for a particular form of 
the perturbation function <p(x). Of course, the same approach could be extended 
to more complicated perturbation functions of the form <f>(x) = 1 — ep(x) for some 
function p(x). The key point consists of determining how the operator correspond- 
ing to the multiplication by p(x) acts on the basic functions hj(x) for j > (defined 
by Equation (jTTJ)). For computing explicit expressions, however, the main difficulty 
is in solving the differential equations satisfied by the coefficients of the expansion. 
When p(x) is a polynomial, a particular solution to the equations similar to Equa- 
tions (|41|) and (|55[) is obtained in the form of a polynomial times the function 
exp(— ot(x — m) 2 /a 2 ) and in that case, explicit computations can be carried out. 

The perturbation function 4>{x) defined by Equation ^ corresponds to the case 
when unresponsive flows have a peak bit rate s much smaller than the transmission 
capacity of the link. The results of this paper show that the reduced service rate 
approximation yields in some conditions accurate results for the performance of 
elastic flows. 



Appendix A. Proof of Proposition [3] 

To prove Proposition [3j we proceed as for the proof of Proposition O We first 
show that the operator (Hi + Ai) -1 ]/^! is bounded, where W\ is the restriction to 
t\(p) of the operator W. (Note that the operator associated with the multiplication 
by — 1 is not bounded in L 2 (R, n).) For this purpose we use the fact that an 
element of / = (fj(x)) G L 2 (R, n; £\ (p)) can be decomposed as 



oo oo 



j=l fe=0 



and the squared norm is 



oo oo 

2 J 



j=l fc=0 

By using the above decomposition, we have 



oo oo 



Bf = / A -MCj,fc + \c j+ltk )h k (x)ej 

j=l fc=0 

By using the recurrence relation satisfied by Hermite polynomials 
(59) xH k {x) = -H k+ i{x) + kH k -i(x), 

we deduce that 



2>N 
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and then 



Wif = e EE G i,khh{x)e } ■, 

3=1 k=0 



where 

Cj,k = m (jtic,-,fc - \c. ]+ i, k ) + 

From the above relations, we have 
(60) {H 1 +A 1 )~ l W 1 f = 



l) + \l —TT- (VCj,k+l - \Cj+l,k+l) 



e J2 C 3,o( H i + A 1 )- 1 h (x)e j + E E c iA H i + A^h^t 

\j=l 3=1 k=l 

For the first term in the right hand side of Equation ([BP]) , we have 



Y,Cj fi {H 1 +A 1 )- 1 h Q {x)e 

3=1 



< Um+A^- 



\(H: 



3=1 

a 



11/11, 



since 



3=1 



\ 



E(<" c i>° - Ac j+ i, )V + --)=, 

i\l OL 



3 = 1 



E^'* 1 - Acj +M )V 

3=1 



together with the inequalities 



EO^.o - Ac,- + i )0 )V < M 



\ 



3=1 



\ 



oo 
3=1 



<m(i + v^)||/|| 



and 



» " Ac 3+i,i)V < m(i + 

For the second term in the right hand side of Equation , since the operator 
H is invertible on the space span(/ifc, k > 1 ), we have 

oo oo oo oo 

E fe ( ffl + M^hkWej = ^ ^ C Jlfc (I + H^ 1 A 1 )- 1 H^ 1 h k (x)e J 

j=l k=l j=l k=l 

and then, by using the fact that ||(I + H^ 1 Ai)~ 1 \\ < 1, we obtain 



J2J2 C ^(H 1 +A 1 )- 1 h k ( 

3 = 1 k=l 



x)e,, 



< 



oo oo / ^ \ 2 



\ ^ E 1 fe 

\ j=i fe=i \ 
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From the inequality 



\ 



oo oo , „ x 2 

EE(%) <* 

j=l k=l V ' 



< m. 



\ 



EX 

3 = 1 k=l 



T 



k 2 



2J^\ EE 

V \ 3=1 fc=l 



k 



2Ja 



(fe + 1) Ocj,fc+i - Ac i+ i, fc+ i) 2 

z^z^ p ^ 



3 = 1 fc=l 



we deduce that 



\ 



OO OO 



J'=l fc=l 



EE(%) ^<M(1 + Vp)(m + ^)||/| 



Finally, by using the fact that + Ai) 1 || < l/(/x(l — v 7 /') 2 ); we come up 

with the conclusion that 



(61) 



(1-^/p) 2 

The operator (ifi + AjJ^Wi is hence bounded. Under Condition the norm 
of this operator is less than 1 and we can adapt word by word the proof of Propo- 
sition^ 
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